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1 Introduction 



Standard cosmology based on the general theory of relativity uses the Friedmann-Robertson- Walker 
(FRW) models, which are perfectly isotropic and homogeneous, or in short 'uniform'. But many 
workers have suggested that the early universe may not have been exactly uniform, especially 
since it must have harbored the inhomogeneities which later evolved to give the galaxies and other 
structures we observe. 

The simplest extension of the FRW models is the anisotropic ones of Bianchi type I, notably that 
of Kasner pQ. In such models, it is possible for the three spatial dimensions to undergo repeated 
"bounces", providing a description for the big-bang fireball if it was initially chaotic and evolved 
towards uniformity. (The use of a 4D Kasner-type metric in this manner leads to the so-called 
mixmaster universe as discussed in ref. 2, but we defer a discussion of the 5D analog to another 
place.) Many anisotropic models for the early universe have been much studied within the context 
of the usual four-dimensional theory of gravity due to Einstein. However, extensions of that theory 
to five and more dimensions are currently seen as providing the best route to unification of gravity 
with the interactions of particle physics [3] . The two versions of 5D relativity in vogue are membrane 
theory and space-time-matter (or induced matter) theory. Both employ a 5D Kaluza-Klein type of 
metric, where the extra dimension is not assumed to be compactified as in the original account and 
is indeed algebraically general [I], [5|. These two modern approaches to 5D relativity are actually 
equivalent in a mathematical sense, but are different as regards physical motivation [BJ . Membrane 
theory allows gravity to propagate in all 5-dimensions, whereas the interactions of particles are 
confined to a singular AD hypersurface, thereby accounting for their relatively greater strength [TJ. 
Space-time-matter theory views the 5D world as pure geometry, where the matter we see in 4D is 
a consequence of the extra dimension in accordance with Campbell's embedding theorem [8], [9]. 
Both theories are in agreement with observations, and both have large literatures of which reviews 
are available [10], [11]. Our present work is applicable to either approach, and is specific: we will 
derive an exact class of solutions to the 5D field equations which describe anisotropic cosmological 
models of Bianchi type-I that may be applicable to the early universe. 

Our solutions will prove to be generalizations of other 5D ones [12]- [IS] and reduce to the 4D 
Bianchi type-I models in the appropriate limits [TB]. The solutions derived bellow generalize, most 
notably, the 'standard' 5D cosmologies [13]. Those involve the fifth coordinate in a significant 
way, but have 4D subspaces with FRW-like properties. They were discovered using the method 
of separation of variables, which ensures that they reduce to the appropriate AD solutions on 
hypersurfaces defined by fixing the fifth coordinate so as to recover spacetime. These solutions 
have one assignable parameter, which determines both the 3D expansion rate and the properties 
of the matter as defined by an effective energy-momentum tensor that describes various forms of 
matter, including radiation (early universe) and dust (late universe). The standard 5D solutions 
also include the analog of the de Sitter solution, which via the cosmological constant describes a 
vacuum-dominated phase (very early universe). However, we will not be much concerned with this 
phase in what follows, because it has recently been the subject of intense investigation [T7]. Instead, 
we will concentrate on the early universe, and use the method of separation of variables to obtain 
exact solutions to the field equations in 5D that have physically acceptable properties of matter in 
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AD. 

For the field equations, we will follow the usual approach and employ the 5D Ricci tensor in the 
form Rab = (A, B = 0, 1, 2, 3, 4). These apparently empty solutions in 5D may contain matter in 
AD, because they include as a subset the Einstein field equations G a p = %nT a p (a,/3 = 0, 1,2,3). 
This is a consequence of the aforementioned theorem of Campbell, which provides a local embedding 
via the formal rules of differential geometry that can be interpreted in a physical way to connect the 
Einstein tensor (Gap) and an effective energy-momentum-tensor (T a p). This method for ensuring 
compliance with the field equations of general relativity has been widely applied. We will employ it 
briefly bellow, but most of our work will be concerned with the "raw" solutions and their associated 
metrics. The latter will be defined via a 5D line element dS 2 = r yABdx A dx B , which contains the 
conventional AD counterpart ds 2 = g a pdx a dx^ . The spacetime coordinates will be assigned as usual, 
with x° = t for time and x 1 ' 2,3 = x,y,z for space. We will write x 4 = ip for the extra coordinate, 
in order to avoid confusion with the spatial Cartesian label and to avoid the implication that x 4 is 
measured with respect to some special hypersurface. The gravitational constant G and the speed 
of light c will be set unity by an appropriate choice of units. Other aspects of our terminology are 
standard. Most of our results are algebraic in nature and derived in Section 2, but we mention some 
astrophysical implications [18j-[21] in Section 3. For rushed readers, equation fflBl defines our class 
of solutions. 

2 Five-dimensional Bianchi Type-I Solutions 

In this section we wish to solve the field equations Rab — using the method of separation of 
variables in x° = t and x 4 = ip. 

We start with a 5D line element where the components of the metric tensor are independent of 
the coordinates for ordinary 3D space (x,y,z) and split via two functions Qa(iP) and SU(i). Thus, 
we write 

dS 2 = Q 2 (^)S 2 (t)dt 2 - Qj(^)S 2 (t)dx 2 - Q 2 2 (^)Si(t)dy 2 - Q 2 3 WS 2 3 (t)dz 2 + eQl{^)Sl{t)d^ 2 . (1) 

Here e = ±1 allows for a spacelike or timelike extra dimension, both of which are physically 
admissible. (See, e.g., references [5], [9], [22]: in modern noncompactified 5D theory, the extra 
coordinate is different in nature from the spacetime coordinates, so there is no problem with closed 
timelike paths). 

Our starting metric (1) is, of course, somewhat special in terms of 5D algebra; but it will be seen 
to be adequate in terms of yielding the appropriate 4D physics. Bearing in mind the opportunity 
to simplify (1), we get rid of the gauge functions S (t) and Q^ip) by means of the coordinate 
transformations 

So(t)dt — > di, and Q^^dtp — > d^. (2) 
Then relabeling the coordinates (dropping the overbars) the metric ([!]) becomes 
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dS 2 = Q 2 {^)dt 2 - Q\{ijj)S 2 {t)dx 2 - Q 2 2 ^)S 2 2 {t)dy 2 - Q 2 {^)S 2 {t)dz 2 + eS 2 4 {t)d^ 2 . 



(3) 



This metric reduces to homogeneous and isotropic Bianchi type-I cosmological models with flat 
spatial sections on every hypersurface : ip = ip = constant. There are eight unknown functions 
in Q, which are to be determined from Rab = 0. 

The field equations have in fact only six nonzero components, which it is convenient to consider 
in the order R 00 , i? 44 , i? 04 , Rn, R22, R 3 3- Setting these to zero gives the following set of equations: 



eSl 



Si So S3 Sa 
—+—+—+— 

Cf Cf Q Q 

Ol 02 ^3 ^4 



+ Ql 



Q'o _|_ Qq_ ( Q\_ _|_ Qj2 _|_ Qji 
Qo Qo \Qi Q2 Q3. 



(4) 



eSl 



Sa Sa I Si . $2 , S3 

q o \ o a a 
04 04 \Ol 02 03 



+ Ql 



Qo 1 Qj_ 1 Q2 1 
Qo Qi Q2 Qs 



(5) 



Q'o I Si s 2 S3 

Qo \Si £2 S$ 



Qi I Si Sa 
Ql \Si Sa 




SA Q'3 ( S3 
Sj Q 3 [S 3 



(6) 



eSl 



Si Si ( S 2 

Q Cf \ Cf Cf 

Oi Oi \0 2 03 




+ Ql 



Qi 1 Qj_ ( Qo_ 1 Qj2 1 Q3 
Qi Qi \Qo Q2 Qs, 



(7) 



eSl 



S2 S2 ( Si S3 S4 
—+——+—+— 

C C 1 Q C C 

O2 O2 Wl ^3 



+ Ql 



Qj 1 Q2 ( Q^o 1 Qi 1 Q3 
Q2 Q2 \Qo Qi Q 3/ 



(8) 



eSl 



S3 S3 ( S\ S2 S4 

Cf ' Q I Q ' Q ' Q 

°3 °3 \°1 °2 04, 



+ Q, 



Q'i _j_ 0^3 ( Qo_ _|_ Qj_ _|_ Q^ 
Q3 Q3 \Qo Qi Q2 , 



0. 



(9) 



In these, a dot and prime mean the total derivative with respect to t and ip, respectively. It should 
be noted that the last three equations are invariant under the interchanges l<-^2, 1^3,2^3, 
as expected. 

We now proceed to construct the most general solution to equations (jl])-©. Towards this end, 
it is helpful to notice a property shared by all of them except the cross-component Rqa given by 
([6]). Namely, that the functions in square brackets which multiply S 2 and Qq must be proportional 
to each other, since otherwise we obtain the trivial result that S4 = and Q = 0. Thus, from ([7]) 
and (IHD we find 
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where k is an arbitrary separation constant which is physically dimensionless (i.e., it is a pure 
number). Further, using the invariance under 1 <-> 2 mentioned before, we obtain from (flQj) the 
relation 



£>2 &2 ( Si S3 S4 
ry ~ Q I O O O 

This may be combined with (fTOj) to give 



rt ^ O I O O O 

Ol Ol \ 02 03 04 



S2 S2 I Si S3 S4 

q ' Q \ Q Q Q 

02 £>2 \ o\ O3 O4 



;i-A; 2 )=0. 



(11) 



(12) 



Since k is an arbitrary number, it follows that 
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(13) 



This relation has analogs that follow from the preceding argument applied to other functions mul- 
tiplying Si, so there are similar functional forms for all the scale functions Si (i = 1,2,3,4). 
Furthermore, (Tl3|) is an ordinary differential equation in x° = t, with no appearance of x A = ip. 
This means that the solution to the field equations takes the usual Kasner form: 

dS 2 = dt 2 - Ai t 2ai dx 2 - A 2 t 2a2 dy 2 - A 3 t 2a3 dz 2 + eA 4 t 2a ^d^ 2 . (14) 

Here the Ai are arbitrary constants whose physical dimensions depend on those assigned to the co- 
ordinates. The ai are dimensionless parameters whose values are constrained by the field equations 
to satisfy the conditions: 



a 1 + a 2 + a 3 + a 4 = 1, 



2222 
a l + a 2 + a 3 + a 4 



1. 



(15) 



Relations like this and its associated metric ([141) also follow when we apply the preceding analytic 
argument not to S"| but to Q^. The analog of the t-dependent solution (FI4|) is the -^-dependent one 



dS 2 = B ^dt 2 - B 1 ij 2l3l dx 2 - B 2 ij 2l32 dy 2 - B 3 ^dz 2 + e d^ 2 
and the analogs of the conditions (TT51) are 

A) + Pi + P2 + /3 3 = 1, K+ >Y+ )'7-h >'; = 1. 



(16) 



/?o + + = 1. (17) 

At this point we have obtained two special Kasner-type solutions, namely the t-form (T141) and the 
■^-form (fl6|) . It remains to combine these to obtain the general solution which we seek. 



Such a general solution of the field equations 
form: 
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Here the constants A — E are arbitrary. The parameters cti, 2,3,4 and /3o,i,2,3 satisfy (1T51) and (ITT)) as 
before; but in addition the cross component i? 4 — of © now imposes an additional condition: 

+ a 2 (3 2 + a 3 p 3 - (3 (ax + a 2 + a 3 ) - a 4 (Pi + P 2 + P3) = 0. (19) 

This condition, together with (IT51) and (ITT)) , might repay a detailed investigation, in order to explore 
the parameter space that can be occupied by the dimensionless constants which define the general 
solution (fl8l) . The latter extends our understanding of Bianchi type-I solutions of Kasner type, in 
the sense of adding to what is known in 4D [JJ, [2], [16] and generalizing what is known in 5D 
[T_2]-[T5]. This is particularly true insofar as ( fl8l) 'overlaps' in form the metric for the standard 5D 
cosmologies [13J. We proceed to make some comments on this and related issues. 

The class of solutions ffl8l) is Ricci-flat (Rab — 0), but in general it is not Riemann-flat (i.e., 
Rabcd 7^ 0). The distinction is important because certain solutions of 5D relativity with high 
degrees of symmetry may have Rabcd = and be flat in 5D, while possessing curved subspaces in 
4D. This is the case for the standard 5D cosmologies [13] . to which we will return shortly. 

In general we note that for our solution ( 1181) a typical component of the Riemann tensor in 5D 
is proportional to (^Mt 2 ijj 2t30 + Nt 2a4 il) 2 ^ , where M and N are some combinations of the parameters 
a and (3. Thus, in general Rabcd 7^ 0. However, for a 4 = (3$ = 1, one can choose the parameters 
in such a way as to make Rabcd = 0. 

Clearly, from (fl5l) and (ITTj) the conditions 0:4 = /3o = 1 imply that a% — a 2 — 0:3 = (3\ — (3 2 — 
/?3 = (assuming that the metric functions are real). However, going back to the field equations we 
find that R\ = R 2 = R^ = requires isotropic expansion for the 3D sections of the metric, which 
entails 

a 1 = a 2 = a 3 = -- fa = fa = Pa = P, e(^)+Ap 2 = 0. (20) 

These values with the Rqa component of the field equation ([6]) imply P = 1/(1 — a), where the 
constant a is defined in (1201) . from which we also get E = —eAa 2 /(l — a) 2 . Drawing together these 
results allows us to rewrite the metric f|T8l) as 

4 2j.2 

dS 2 = Ai?dt 2 - B^' a k 2 / a (dx 2 + dy 2 + dz 2 ) - -j^^dip 2 . (21) 

v ' (1 — a) z 

It is important to note that a is a free parameter, i.e., it is not constrained by the field equations. 
The constants A and B can in principle be absorbed by a suitable choice of units, when the metric 
becomes identical to that for the standard 5D cosmologies with flat space sections [13]. These models 
have been much discussed (see [3], p. 35 for a review). On the hypersurfaces x 4 = ip = constant, 
they give back the conventional FRW models for the early radiation- like universe (a = 2) and the 
late dust-like universe (a = 3/2). These properties of matter, and other aspects of the solution 
f l2T|) . are very special in the present context, however, because that case is the only member of the 
general class fTTSl) which is Riemann-flat. 

The properties of matter for the general metric ffTSl) may be evaluated using the standard 
technique based on Campbell's theorem [S], [9]. In the present situation this consists in isolating 
the 4D part of the relevant 5D geometric quantities, using them to construct the 4D Einstein tensor 
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G af 3 (a,/? = 0,1,2,3) and using the field equations of general relativity G a/ 3 = 8irT al 3 to identify 
the effective 4D energy-momentum tensor (see section 1). The working is tedious, so we just quote 
the results: 



87rT ( ° 
8ttTI 
87t71 



oi(o 2 + o 3 ) + a 2 a 3 
A^H 2 ' 

a 2 + a l + «2(«3 - 1) — a 3 

2 of + «3 + o x (o 3 — 1) — o 3 
Aip 2 P°t 2 ' 



g7rT 3 = a; + al + ai(a 2 -l)-a 2 ( } 

Here the parameters are constrained by (}T5]) . (fT7|) and f|T9|) . We note that the last three components 
are invariant under the interchanges l<->2,1^3,2^3,a property that they inherit from the field 
equations (jTJ), (jSJ) and (jUJ). Another property of fl22l is that, for a general choice of parameters, the 
components satisfy T = T$ + T 1 + Tf -t-Tf = 0. This in 4D physics is usually taken to mean that the 
equation of state is that of radiation (i.e., photons with zero rest mass) or ultra-relativistic matter 
(i.e., particles with finite rest mass moving close to the speed of light). However, for the so-called 
induced-matter theory, it is well-known that when the 4D part of the 5D metric is independent 
of x 4 , the effective energy-momentum tensor has T = (see e.g. [3]). Our result implies that the 
reverse implication need not hold. We also note that the isotropic cosmology (1201 has 



T = 2(2-a)T °, 87rT ° = -^A_ T { = Tf = Tf = (l - |) T °. (23) 

Since a is an arbitrary parameter we have T^O, except for radiation (a = 2), as expected. 

From fl22l it follows that T\ = T 2 = T 3 3 requires either 0:1 = 0:2 = o 3 or 01 + o 2 + o 3 = 1. But 
the former choice leads to the isotropic solution (j2"Tp . while the latter choice implies o 4 = by ffTol) 
and leads back to the 4D vacuum Kasner solution pQ, [2]. This means that the 4D effective matter 
source for the 5D solution ([TBI cannot in general be modeled as a perfect fluid. This is consistent 
with other comments on the comparable 4D solutions. For example, one authoritative text refers 
to the related mixmaster cosmology as showing "properties of empty space reminiscent of an elastic 
solid" (see [2j p. 806). 

In addition to providing 5D analogs of the 4D Kasner solution, our metrics may also be related to 
others in the braneworld scenario [3], [1], [7], [10]. In that approach there is a singular hypersurface 
that defines spacetime, and the properties of matter in that hypersurface may not be identical to the 
effective ones we have calculated above from the induced-matter approach as based on Campbells 
theorem (8] , [9] . Our solutions may be verified by the application of a program such as GR- Tensor 
to satisfy Einsteins equations. However, it is well known that in higher-dimensional physics a 5D 
metric may under coordinate transformations correspond to different 4D sources. For this and 
other reasons, we recommend further investigation of the matter and dynamics associated with our 
general class (ITHI) of solutions. 
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3 Summary 



We have derived exact solutions of the five- dimensional field equations which extend those of general 
relativity that are of Bianchi type-I. Our general solution ( fl8l) has a number of parameters that are 
constrained by relations (fl5l) . (TlTI) and ([191) . It has similarities with the Kasner metric of Einstein's 
theory, and can be regarded as a higher-dimensional version of that model. 

It is logical to apply our solutions to the era after the big-bang, because the solutions (TTSl 
have a 4D source (1221) with a radiation-like equation of state, and because the special form (12T1) is 
known to provide a reasonable description of the universe at later epochs [3], [ 1 3J . We could, for 
example, imagine a model which proceeds through a series of physical phase changes, the phases 
corresponding to different values for the assignable parameters in the metric ffl8|) . One such model 
has already been investigated [18], in order to explain the transition from deceleration to acceleration 
suggested by recent supernova data jT9] - Other consequences of 5D models have been studied [20J, 
in order to avoid some of the problems to do with the formation of galaxies and other structures in 
conventional 4D cosmology [2T] . We have not, for want of space, investigated these issues here, but 
they are worthy of attention. Likewise, it has been suggested that particles might be produced in 
the early stages of an anisotropic universe [2] , and the scalar field represented by the last component 
of a 5-D metric would provide a natural mechanism for matter production. 
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